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ABSTRACT 

We present full stellar evolution and parametric models of the surface abun- 
dance evolution of 16 0, 22 Ne, 23 Na, and the magnesium isotopes in an extremely 
metal-poor intermediate mass star (A^zams = 5 Z = 0.0001). 16 and 22 Ne 
are injected into the envelope by the third dredge-up following thermal pulses on 
the asymptotic giant branch. These species and the initially present 24 Mg are 
depleted by hot bottom burning (HBB) during the interpulse phase. As a result, 
23 Na, 25 Mg and 26 Mg are enhanced. If the HBB temperatures are sufficiently 
high for this process to deplete oxygen efficiently, 23 Na is first produced and then 
depleted during the interpulse phase. Although the simultaneous depletion of 
16 and enhancement of 23 Na is possible, the required fine tuning of the dredge- 
up and HBB casts some doubt on the robustness of this process as the origin 
of the O-Na anti-correlation observed in globular cluster stars. However, a very 
robust prediction of our models are low 24 Mg/ 25 Mg and 24 Mg/ 26 Mg ratios when- 
ever significant 16 depletion can be achieved. This seems to be in stark contrast 
with recent observations of the magnesium isotopic ratios in the globular cluster 
NGC 6752. 

Subject headings: stars: AGB — stars: abundances — globular clusters: general 



1. Introduction 



In globular clusters (GCs) spectroscopic observations have revealed large (~ 1 dex) star- 
to-star abundance variations of C, N, O, Na, Mg and Al (e.g. Ramirez & Cohen 2002). The 
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anti-correlations of C-N, O-Na and Mg-Al point to simultaneous operation of the CNO-, 
NeNa-, and MgAl-cycles. In the so-called evolutionary scenario it is assumed that these 
abundance variations are produced in the vicinity of the hydrogen-burning shell in red giant 
branch (RGB) stars and that some extra-mixing transports them to the convective envelope 
(Denisenkov & Denisenkova 1990; Langer et al. 1993; Denissenkov & Weiss 2001). 

The evolutionary scenario has been challenged by the recent discovery of C-N, O-Na 
and even Mg-Al anti-correlations in the main sequence (MS), MS turn-off and subgiant stars 
in the GCs 47 Tucanae and NGC 6752 (Harbeck et al. 2003; Gratton et al. 2001; Grundahl 
et al. 2002). A likely origin of these abundance anomalies is pollution by material processed 
via H-burning in more evolved stars. In this primordial scenario intermediate-mass stars 
(IMS, M. « 3-8.M©) in their asymptotic giant branch (AGB) evolution phase have been 
proposed as contaminators (Ventura et al. 2001, 2002). 

Thermally pulsing (TP-) AGB stars process material in a He- and a H-burning shell. 
The H-shell dominates the energy production most of the time. However, recurrent ther- 
monuclear He-shell flashes drive a temporary (~ 10 yr) pulse-driven convective zone (PDCZ) 
that encompasses the entire region between the He- and the H-burning shell (the intershell). 
Immediately after the end of the TP the base of the convective envelope begins to move in- 
wards in mass and eventually material from below the H-shell is dredged-up to the envelope 
(third dredge-up, TDU). During the interpulse period in massive AGB stars, hot-bottom 
burning (HBB) further modifies the envelope chemical composition (Boothroyd et al. 1993). 

In the evolutionary scenario Al can be synthesized on the RGB if some low-mass stars 
in GCs were initially enriched with 25 Mg, possibly from massive AGB stars (Denissenkov 
et al. 1998; Denissenkov & Weiss 2001). In this combined scenario star-to-star abundance 
variations of C, N, O, Na, Mg and Al in GCs may have multiple origins: (i) on the MS they 
may be due to pollution in the past either by the massive AGB stars or by somewhat more 
massive (0.9 < M.j 'A4® < 2) RGB stars than the present-day MS turn-off stars which had 
undergone deep extra- mixing, and (it) on the RGB they may be partly tracers of the same 
pollution that occured on the MS and partly (in the most rapidly rotating stars) they may 
be due to deep extra-mixing (Denissenkov & Weiss 2001; Denissenkov & VandenBerg 2003). 

In the primordial scenario massive AGB stars are thought to be responsible for very low 
O abundances ([O/Fe] < —0.5, down from the assumed initial value of +0.4) in MS stars in 
GCs by pollution with O-depleted material. Indeed, Ventura et al. (2001, 2002) have reported 
that, in metal-poor massive AGB stars, HBB may be capable of producing the required O 
depletion. However, nucleosynthesis beyond the CNO-cycle has not been investigated yet. 
In this Letter we will take into account all of the cycles of nuclear reactions participating in 
the H-burning as well as the effect of the TDU. In § 2 we present the abundances from full 
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stellar evolution models. In § 3 an equivalent parametric AGB model is described. In § 4 we 
discuss the results of calculations with the parametric models and make our final conclusions. 

2. Full Stellar Evolution Models 

Our ID stellar evolution code (Herwig 2000) includes updated opacities (Iglesias & 
Rogers 1996; Alexander & Ferguson 1994) and a nuclear network with all relevant reactions 
with the rates for the NeNa-cycle reactions taken from El Eid & Champagne (1995). Simul- 
taneous, fully implicit, iterative solution of the nuclear network and time-dependent mixing 
equations for each isotope (Herwig 2001), and hydrodynamic overshooting with a geometric, 
exponential decay parameter / can be included. We have improved the adaptive time-step 
and grid-allocation algorithm for the extremely metal-poor models in order to ensure that 
the TDU properties are not affected by numerical resolution issues. The mixing-length pa- 
rameter is «mlt = 1-7 from calibrating a solar model, and X( 4 He) init = 0.23. 

We choose a model of initially 5Ai & with a metallicity of Z = 0.0001 as a representative 
example for IMS that might have polluted GC stars of the lowest metallicity ([Fe/H]f« 
\og(Z/Z & ) = —2.3). We evolve the initial pre-MS model through all evolutionary phases, and 
details will be presented in a forthcoming paper. Up to the first thermal pulse on the AGB, 
exponential overshooting with / = 0.016 has been considered at all convective boundaries. 
The first envelope abundance alteration occurs as a result of the second dredge-up. Material 
processed mainly by H-shell burning is brought to the surface. 23 Na is enhanced by 0.75 dex 
from the conversion of the initial abundance of 22 Ne and some 20 Ne. 16 is depleted by 
0.1 dex due to the action of the ON cycle. Magnesium isotopes are changed by less than 
0.05, 0.1 and 0.01 dex for mass numbers 24, 25 and 26 in the second dredge-up. 

During the TP-AGB phase the interplay of the TDU and HBB alters the envelope 
abundances. TDU brings 22 Ne and a rather uncertain amount of 16 to the surface. 25 Mg 
and 26 Mg are produced by two processes: (i) by HBB in the envelope (at the expense 
of 24 Mg), and (ii) by a captures on 22 Ne in the PDCZ. The second process is important 
when the temperature in the PDCZ exceeds ~3.5 x 10 s K. The HBB also destroys 16 (and 
produces 17 0), while 23 Na is produced by proton captures on the dredged- up 22 Ne. At higher 
temperatures, 23 Na can be destroyed again. The overall budget of the O, Na and Mg isotopes 
depends on the interplay of dredge- up and HBB. In addition, mass loss plays an important 
role. It limits the TP-AGB evolution time with high HBB temperature and efficient third 
dredge- up. 

All of these effects can be observed in the surface abundance evolution of two TP-AGB 
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model sequences that we evolved from the same early- AGB model (Fig. 1). In the first 
sequence (a) we have isolated the effect of HBB by assuming no overshooting and a rather 
low time resolution (typically less than 1000 models per TP cycle). As a result, this sequence 
shows no TDU (apart from one outlier seen in the run of 16 at t — 60000 yr). We have also 
assumed mass loss according to Blocker (1995) with rjs = 0.1, and the model experiences 43 
TPs until the envelope mass is lost. This sequence shows a large 16 depletion, but due to 
the lack of 22 Ne dredge-up, 23 Na is depleted as well. 

In panel (b) a model is shown in which we have considered only negligible overshooting 
at the bottom of the He-flash convection zone (/ = 0.002; the intershell abundances are not 
affected significantly by this very small PDCZ overshoot). At the bottom of the envelope 
convection we have assumed / = 0.016, as calibrated at the MS core convection boundary. 
Here we have used the high spatial and time resolution required to model the TDU. Mass loss 
has been turned off for this case in order to explore the largest possible nuclear enrichment 
by HBB and dredge-up. After a few TP the TDU is very efficient, and it reaches into 
16 0-rich layers below the zones previously covered by the PDCZ. Our models with this deep 
dredge-up show systematically higher HBB temperatures than the model without dredge-up. 
The unusually deep dredge-up is likely related to the effect of sustained, fierce H-burning in 
the overshoot layer. We will discuss the details of this effect, which is unique to very metal 
poor IMS, in a forthcoming paper. For the present study it is sufficient to note that the 
TDU in these stars could be much more efficient than previous models without overshooting 
indicated (e.g. Ventura et al. 2002), and that oxygen could be enhanced by this process. This 
undermines the ability of HBB to deplete oxygen in the envelope. However, with efficient 
TDU, 22 Ne is dredged-up and 23 Na is, in fact, further increased. The behaviour of the 
magnesium isotopes is qualitatively the same as in the model without dredge-up, indicating 
that these species are mainly affected by the HBB. 

According to our high-resolution, full stellar evolution models with efficient TDU, mas- 
sive AGB stars cannot show simultaneously 16 depletion and Na enhancement. Using 
parametric models we will now explore intermediate cases with 22 Ne dredge-up but limited 
16 dredge-up, and we will consider in more detail the evolution of magnesium isotopes. 

3. Parametric Models 

We have developed a parametric nucleosynthesis and mixing code for massive AGB 
stars. All processes that are relevant for the abundance evolution of O, Na, and the Mg 
isotopes are considered. In particular, we include the effects of HBB and dredge-up, as 
well as mixing and burning in the PDCZ. Most reaction rates are taken from the NACRE 
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(Angulo et al. 1999) compilation (see Denissenkov & Tout 2003, for details). For the initial 
conditions we assume scaled solar abundances except for an enhancement of the a-elements 
( 16 0, 20 Ne, 24 Mg, etc.) by +0.4 dex, and a depletion of Na and Al by the same factor, in 
accordance with the chemical composition of the halo dwarfs. As in the full stellar models, 
the metallicity is Z = 0.0001. In the following description, we use the subscripts "PDCZ", 
"HBS" and "HBB" to indicate the temperatures and densities at the base of the PDCZ, in 
the H-burning shell and at the bottom of the convective envelope, respectively. 

A TP sequence starts with a calculation of the mixing and nucleosynthesis in the PDCZ. 
The pre-TP composition of the intershell zone is a mixture consisting of a fraction 7 of 
material from the preceding PDCZ and a fraction (1 — 7) containing H-shell ashes. For the 
first TP we take 7 = 0. The mass AA^pdcz of the PDCZ is divided into 25 zones with the 
temperature and density decreasing linearly from Tpdcz and Ppdcz at its bottom to Thbs 
and phbs at its top. For the mixing in the PDCZ, which is treated as a diffusion process, we 
assume a constant coefficient D mix = 10 15 cm 2 s" 1 . For the temperature T PD cz we choose the 
maximum value observed in the full stellar models. For that reason we are really modeling 
only the high-temperature phase of the PDCZ. The He-burning is stopped every time when 
the mass fraction abundance of 12 C in the PDCZ has reached the value 0.23, as proposed 
by Renzini & Voli (1981). After that a fraction A of the total mass AA^pdcz of material 
with the final abundances from the PDCZ is added to the envelope to simulate the effect of 
the TDU. Then, we follow the changes of the envelope composition due to the HBB during 
the interpulse period At ip . For this, we process the envelope abundance distribution after 
the TDU in one zone at constant T H bb and Phbb<5 (a dimensionless factor S accounts for 
the fact that, after averaging over the mass of the convective envelope, coefficients in the 
nuclear kinetics network can be written in a form (av )(Thbb)-^aPhbb<5, where (av)(T)N\ is a 
reaction rate at a temperature T and 5 = (1 — a) -1 j a (p/ p-rbb)[{<7v) (T) / (cry) (Tuw)]dx with 
x = M r /M and a = Mbbb/M; the full stellar evolution models give 5 ~ 10~ 5 , which means 
that HBB takes place in a narrow zone adjacent to the base of the convective envelope). In 
addition to the HBB computation, we calculate the abundance distribution of the H-burning 
ashes from the new post-HBB envelope composition with another one-zone model at H-shell 
temperature and density. This whole sequence is repeated for eight TP cycles. 

We use the following structure parameters from the full models: the CO-core mass 
M c = O.96A4 , AAWz = 0.0025 M Q , T PDCZ = 3.2 x 10 8 K, p PDC z = 5 x 10 3 gcm- 3 , 
Thbs = 1-2 x 10 8 K, p HBS = 40gcm~ 3 , T HBB = 10 8 K, p HBB 5 = lO^gcm" 3 , and At ip = 
10 4 yr. Our test calculations have shown that the envelope abundances depend weakly on 
the parameter 7 for a wide range of 0.2 < 7 < 0.6. Therefore, its value was kept constant 
at 7 = 0.6, as given by the full stellar evolution calculations. 
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The efficiency of nuclear processing of the envelope material in the HBB is determined 
by 3 parameters: At ip , Phbb<5 and T H bb- However, while the efficiency is linearly proportional 
to the first two of them, it depends on a high power of Thbb- Therefore, in our model only 
the third parameter characterizes the efficiency of HBB. The envelope abundances are also 
strongly affected by the efficiency of the TDU (parameter A). In the full stellar evolution 
models A depends critically on the efficiency of convection induced extra-mixing, while Thbb 
depends on the efficiency of convective energy transport in the envelope. In our parametric 
model A and T H bb are free parameters. 



4. Results from the parametric model and general discussion 

Parametric calculations were carried out for A = 0.3 and 1, and Thbb = 0.9 ... 1.1 x 10 8 K. 
The time evolution of the envelope abundances (Fig. lc) resembles the main features seen in 
the results of the full stellar evolution calculations: the dredge-up of 22 Ne, production and 
destruction of Na via 23 Na(p,a) 20 Ne and 23 Na(p,7) 24 Mg, depletion of 24 Mg and production of 
25 Mg and 26 Mg. 16 is depleted just by HBB. In our standard parametric model we assume 
that 2% of the 16 O mass fraction is dredged- up from the CO-core during every TP, in 
agreement with the average minimum predicted by the full stellar models. In the parametric 
model the 16 abundance becomes significantly depleted with Xfi n ( 16 0)/X init ( 16 0) < 0.2 for 
Thbb > 10 8 K and A = 0.3 (Fig. 2). Thus, if no or very little 16 is brought to the surface, 
then we confirm the findings of Ventura et al. (2001, 2002) that 16 is efficiently destroyed 
by the HBB in the metal-poor massive AGB stars. However, our parametric models confirm 
the result of the full models that high HBB temperatures do not favour Na production 
which would be required to explain the O-Na anti-correlation in GC stars in the primordial 
scenario. Indeed, at T = 10 8 K, the sum of the rates of the reactions 23 Na(p,a) 20 Ne and 
23 Na(p,7) 24 Mg is ~6.8 times as large as that of 16 0(p,7) 17 F, the lower and upper limits of 
this ratio being 2.1 and 45, respectively (Angulo et al. 1999). Without the 22 Ne dredge-up 
source which replenishes Na (as in our first full stellar evolution model), the O depletion is 
accompanied by Na destruction at T = 10 8 K. 

But even with the dredge-up of 22 Ne, the final Na abundance in the envelope is found to 
be very sensitive to small variations of Thbb: increasing Thbb from 0.9 x 10 8 K to 1.1 x 10 8 K 
turns the Na production into the Na destruction (Fig. 2). Consequently, for high HBB 
temperatures, both O and Na may be depleted simultaneously, resulting in O-Na correlation 
instead of the O-Na anti-correlation. If dredge-up of 16 O is as efficient as predicted by the 
high-resolution full stellar models with hydrodynamic overshooting, the temperatures for 
efficient 16 depletion would definitively be too high for Na production. Furthermore, Iliadis 
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et al. (2001) have recommended to use a rate of the reaction 22 Ne(p,7) 23 Na which is by a 
factor of ~ 10 2 - 7 smaller than the value given by NACRE at T = 10 8 K. This makes Na 
production in massive AGB stars even more problematic. On the other hand, a rate of 
the reaction 22 Ne(p,7) 23 Na from El Eid & Champagne (1995), used in §2, is by a factor of 
~ 10 3 ' 9 larger than the NACRE value, which favours Na production. In summary, we cannot 
entirely rule out that contamination by massive AGB stars causes the O-Na anti-correlation 
in the primordial scenario. However, this seems to be very unlikely because it requires a 
fine tuning of the AGB model parameters. This is not supported by the latest rate of the 
reaction 22 Ne(p,7) 23 Na either. 

A very robust prediction of the primordial scenario can be made with respect to the 
magnesium isotopes. For temperatures in the vicinity of T H bb ~ 10 s K, that allow 16 deple- 
tion, the 24 Mg abundance is depleted in any case even more, producing 25 Mg and 26 Mg in 
turn. This result is consistently found in all of our models for a wide range of parameters. 
The evolution of 24 Mg in relation to 16 is a robust result because: (i) at T = 10 8 K, the ratio 
of the 24 Mg(p,7) 25 Al rate to the 16 0(p,7) 17 F rate is ~ 6.4, the lower and upper limits being 
3.9 and 13 (Angulo et al. 1999), and (ii) there is no source of 24 Mg in AGB stars (production 
of 24 Mg in the reaction 23 Na(p,7) 24 Mg in the HBB is unimportant). 

Therefore, in the primordial scenario with massive AGB stars as the contaminators, 
the MS turn-off stars with [O/Fe] < —0.5 must have 24 Mg depleted and 25 Mg enhanced 
by more than one order of magnitude. Our assumed initial chemical composition has 
[O/Fe]=0.4, [ 24 Mg/Fe] = 0.4, [ 25 Mg/Fe] = 0.0 and [ 26 Mg/Fe] = 0.0, i.e. the Mg isotopic 
ratios 24 Mg : 25 Mg : 26 Mg = 90.5:4.5:5.0. According to those models in Fig. 2 with sig- 
nificant O-depletion, material released by massive AGB stars may have log(X fin /X ini t) « 
-1.0, -1.5, 1.2 and 0.5 for O, 24 Mg, 25 Mg and 26 Mg (in Fig. 2 the label 26 Mg represents the 
sum of 26 Mg and 26 Al). In low-mass RGB stars in the primordial scenario, ~ 90% of the mass 
of the convective envelope has to consist of this material. Accordingly, on the upper RGB, 
such a star would have [O/Fe] = -0.32, [ 24 Mg/Fe] = -0.50, [ 25 Mg/Fe] = 1.2, [ 26 Mg/Fe] = 0.5, 
and the Mg isotopic ratios 24 Mg : 25 Mg : 26 Mg = 13 : 71 : 16. 

These predictions seem to be in conflict with the results of Mg isotopic composition 
analysis of RGB stars in the GC NGC 6752 reported by Yong et al. (2003). In the least 
polluted stars these authors infer 24 Mg : 25 Mg : 26 Mg^ 80 : 10 : 10 and [O/Fe] ^0.6 while in 
the most contaminated stars they find 24 Mg : 25 Mg : 26 Mg^ 60 : 10 : 30 and [0/Fe]?» —0.1. 
Thus, despite having O depleted by a factor of ~ 5, the NGC 6752 red giants still exhibit the 
24 Mg dominated isotopic ratios. Moreover, the second most abundant isotope is 26 Mg instead 
of 25 Mg. We see only two possibilities to remove this disagreement within the primordial 
scenario: either (i) the ratio of the reaction rates of 24 Mg(p,7) 25 Al and 16 0(p,7) 17 F at 
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T ps 10 8 K is much less than the value given by Angulo et al. (1999), and at the same time 
the reaction 25 Mg(p,7) 26 Al is faster; or (ii) the HBB temperature in massive metal-poor AGB 
stars is somewhat lower than 10 8 K, in which case the 24 Mg destruction would be suppressed 
(see the open circles in Fig. 2). However, in the second case, 16 destruction would be 
suppressed as well, and 25 Mg (and not 26 Mg) could still be produced in a large amount. 
Therefore, we would need deep extra-mixing on the RGB to deplete O and to produce 26 Mg 
and probably Al (either in the form of 27 A1 or 26 Al g ) at the expense of 25 Mg. This is exactly 
the combined scenario proposed by Denissenkov et al. (1998), with the minor modification 
that deep extra-mixing in the RGB stars slightly more massive than the present-day MS 
turn-off stars in GCs might have contributed to the star-to-star abundance variations as 
well. 
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Fig. 1. — Envelope abundance evolution for two full TP-AGB stellar model sequences with 
-Mzams — 5A4 Q , Z = 0.0001 (panels a and b). The sequences are distinguished by differing 
assumptions concerning the mixing and numerical time resolution, as discussed in the text. 
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Fig. 2. — Final envelope abundances with respect to the initial ones after 8 thermal pulses 
calculated with the parametric AGB model: (A, Thbb) = (1, 10 8 K) {filled squares), (0.3,9 x 
10 7 K) {open circles), (0.3,10 8 K) {asterisks), and (0.3,1.1 x 10 s K) {filled circles). 



